Quantum phases of a Feshbach-resonant atomic Bose gas in one dimension 
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We study an atomic Bose gas with an s-wave Feshbach resonance in a one-dimensional optical 
lattice, with the densities of atoms and molecules incommensurate with the lattice. At zero temper- 
ature, most of the parameter region is occupied by a phase in which the superfluid fluctuations of 
atoms and molecules are the predominant ones, due to the phase fluctuations of atoms and molecules 
being locked by a Josephson coupling between them. When the density difference between atoms 
and molecules is commensurate with the lattice, two additional phases may exist: the two compo- 
nent Luttingcr liquid where both the atomic and molecular sectors are gapless, and the inter-channel 
charge density wave where the relative density fluctuations between atoms and molecules are frozen 
at low energy. 
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Introduction Trapped dilute cold atomic gases are 
one of the most exciting fields in condensed matter 
physics. An important recent development in this area 
is the application of Feshbach resonances. The energy 
difference between the molecular state and the two-atom 
continuum, known as the detuning v, can be experimen- 
tally tuned by means of a magnetic field. Therefore, by 
sweeping the magnetic field from positive to negative de- 
tuning through the Feshbach resonance, it is actually pos- 
sible to form molecules in the atomic gas.Q 

Depending on the quantum statistics of atoms, the low 
temperature properties of a dilute atomic gas with an 
s-wave Feshbach resonance will be much different. In 
the case of fermonic atoms, a crossover from a super- 
fluid (SF) of the Bardeen-Cooper-Schrieffer type to the 
Bose-Einstein condensation (BEC) is predicted. 0] On 
the other hand, a quantum phase transition can occur 
for bosonic atoms by changing the value of magnetic 
field detuning. There, two thermodynamically distinct 
phases exist at zero temperature: the "atomic superfluid" 
(ASF) phase with both atomic BEC and molecular BEC 
and the "molecular superfluid" (MSF) with molecular 
BEC only. 
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FIG. 1: Schematic phase diagram of a Feshbach-resonant 
atomic Bose gas in a one-dimensional optical lattice: (a) 
on-site interactions only and (b) Tonks limit with nearest- 
neighbor repulsions between atoms (molecules) V a (V m ). 8 
denotes the density difference between atoms and molecules. 
ao is the lattice spacing. 



In the present letter, we study the possible quantum 
phases for an atomic Bose gas with an s-wave Feshbach 
resonance in a one-dimensional optical lattice. It can be 
described by the Hamiltonian 
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Here n a ^ i is the atom (molecule) density operator and /i 
is the chemical potential. This system can also be consid- 
ered a kind of binary mixtures with two types of bosons. 
However, a main distinction between the present case 
and the binary mixture studied previously |5|,|6( is that a 
Josephson coupling between atoms and molecules, the U 
term in Eq. is absent in the latter. As we shall see, 
this will change the global phase diagram dramatically. 

Our main results can be summarized in Fig. ^ For on- 
site interactions only, i.e. the Hamiltonian given by Eq. 
(|T|). the Josephson coupling dominates the low-energy 
physics in most of the parameter region, and the result- 
ing phase exhibits strong superfluid fluctuations in both 
atomic and molecular sectors. This is because the phase 
fluctuations of both sectors are locked by the Joseph- 
son coupling. This phase is, in fact, a one-dimensional 
(Id) analog of the ASF phase in three dimensions (3d). 
(But we do not have real condensates here.) Therefore, 
it will be referred to as Id ASF. The Id ASF will col- 
lapse for sufficiently strong attractions between atoms 
and molecules (U a b < 0). When the density difference 
between atoms and molecules is close to some ratio- 
nal number, depending on the values of parameters as 
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shown in Fig. ^ (b), two additional phases may emerge 
by including nearest-neighbor repulsions between atoms 
(molecules): the two component Luttinger liquid (2-LL) 
where both the atomic and molecular sectors are gap- 
less, and the inter-channel charge density wave (I-CDW) 
where the relative density fluctuations between atoms 
and molecules are frozen at low energy. 

The continuum theory We now outline below the 
derivation of our results. For simplicity, we shall con- 
sider the case where both the densities of atoms and 
molecules, p a = (n a< i)/ao and p — (n^/ao, respec- 
tively, are incommensurate with the lattice. (Here ao 
is the lattice spacing. p a and pi, satisfy the constraint: 
p a +2pi, = po where po is the density of bare atoms.) That 
is, we do not consider the possibility of the SF-Mott in- 
sulator transition. Then, in terms of the "bosonization" 
formula^: a 4 ~ e -W^^(x) J^L^ e <2™ Pa x e is/ton<t> a (x) 
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low-energy physics of H [Eq. QJ] can be described by 
the following effective Hamiltonian: 
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and Luttinger liquid (LL) parameters, respectively. In 
Eq @, only those terms which may become the most 
relevant in the renormalization group (RG) sense are re- 
tained. The values of v a / m , K a / m , gi, g 2 , and g% depend 
on the short-distance physics. In general, they must be 
extracted from numerics or experiments. K a / m 3> 1 in 
the weak coupling regime. On the other hand, K a / m = 1 
in the Tonks limit, i.e. f7 (-ro) Aa(m) — > +00. Therefore, 
for on-site interactions only, 1 < K a / m < +00. || The 
value of K a (K m ) can be further decreased by including 
nearest- neighbor repulsions between atoms (molecules). 

The g 3 term is a Umklapp process. It can be neglected 
when S is not close to zero (incommensurate filling). On 
the other hand, when 5 is close to zero (commensurate 
filling), the 173 term can affect low-energy physics and 
one may no longer neglect it in Eq. (J2J). When ao|<5| is 
close to some rational number k/l where k = 1,2, 
I = 2, 3, ■ ■ ■, and k and / are co-prime with one another, 
other Umklapp processes must be taken into account, 
and one must include the following term in Eq. J2J: 



We shall focus on Eq. (J2J) in the following and discuss 
the effects of the 53 term later. 

Incommensurate filling When i5 is incommensurate 
with the lattice, one may neglect the 53 term in Eq. 
©. To analyze the effects of the g\ and g 2 terms on 
low-energy physics, we resort to the RG method. A per- 
turbative calculation up to the one-loop order already 
shows that the g\ and g 2 terms alone do not form a 
closed operator algebra in the sense of operator product 
expansion (OPE). One must include the term d x 9 a d x m 
in H e ff [Eq. (J2J]. To simplify the analysis, we con- 
sider the case v a — v m = vq. (The effects of velocity 
anisotropy will be discussed later.) Further, we rescale 
gi by a^gi — > g±. The one- loop RG equations can be 
obtained by setting X 3 (l) = in Eqs. JUJJ — l(T5jl. 
By solving these scaling equations, one may obtain the 
following results: Within the weak-coupling region, i.e. 
ir\gi\/(2v ), \g 2 \/(2vo) -C 1, the g\ term is relevant in the 
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Therefore, there are two zero-temperature phases. It 
turns out that the phase transition between these two 
phases is of the KT type. We note that for on-site inter- 
actions only, i.e. 1 < K a / m < +00, the g\ term is always 
relevant from Eq. 

In the regime where the g\ term is relevant, it is con- 



venient to define new bosonic fields: 6± 
\<i>n ± <t>m.- Then, the value of (9- 
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} is pinned and a 

gap A ~ [7r|(/i |/ (2z; )] 2— jD;L is opened for the 9- sector. 
That is, the phase fluctuations of atoms and molecules 
are locked by the Josephson coupling. By integrating out 
the gapped sector, the low-energy effective Hamiltonian 
describing the gapless (9 + ) sector takes the form of LLs: 
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For 7r|<7i|/ (2vq), |<72|/(2t>o) and K can be related to 

the short distance variables with the help of the one-loop 
RG equations, yielding 
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FIG. 2: Phase diagram for S — 0: (a) suggested by tree- 
level scaling analysis and (b) predicted by the one-loop RG 
equations. Di < 2 and Di > 2 in Region I, Dx > 2 and 
D 2 < 2 in Region II, D 1 ,D 2 >2 'm Region III, and Di, D 2 < 2 
in Region IV. Point A and B in (b) are tricritical points. 



is because v < 0, indicating an instability of this system. 
In other words, the system will collapse for sufficiently 
strong attraction between atoms and molecules. 

The Id ASF can be characterized by the single-particle 
Green functions of atoms and molecules: 



(7) 



where x = (r,x), r = x/(v+t) 2 + x 2 , 9 a (x) = cu/y/ao, 
^m(x) = bi/^/Eo, A\ t 2 are nonuniversal constants, and 



cto = 4ai = 

2K 



We would like to stress that this exact relation between 
ot\ and «2 [Eq. JjjJ] results from the fact that the phase 
fluctuations of atoms and molecules are locked by the 
Josephson coupling, and is a characteristic of the Id ASF. 
On the other hand, the density correlation functions of 
atoms and molecules at nonzero momenta decay expo- 
nentially due to the presence of the gap A . Equation 
|(SJ implies that both the atomic and molecular sectors 
exhibit the behavior of a Id SF as long as K > 1/4. For 
1 < K a i m < +oo, this condition is satisfied, which can 
be verified from Eq. ©. It gives further support to our 
claim that this phase is a Id analog of the ASF in 3d. 

In the regime where the gi term is irrelevant, the ex- 
citation spectrum consists of two branches of gapless ex- 
citations with linear dispersion relations. This phase, 
which has been thoroughly discussed in Ref. 0, will be 
referred to as the 2-LLs. We just mention that in the 
present situation this phase can exist only by including 
sufficiently strong repulsions between atoms or molecules. 
This is because D\ < 2 for on-site interactions only. 

Commensurate filling When S is close to zero, the g 3 
term must be retained. For simplicity, we shall consider 
6 = 0. Let us first define the scaling dimension of the g 3 
term: 



The gi (g 3 ) term will be relevant if Di < 2 (D 2 < 2). Ac- 
cordingly, a tree-level scaling analysis suggests the phase 
diagram in the K a -K m space as shown in Fig. [21 (a). 
D\ < 2 and D 2 > 2 corresponds to the Id ASF, D x > 2 
and D 2 < 2 will be referred to as the I-CDW, and 
D\,D% > 2 corresponds to the 2-LL. However, a com- 
petition between the two relevant operators, the gi and 
33 terms, occur when D\,D2 < 2. To determine whether 
Region IV in Fig. (a) corresponds to a new phase or 
not, we employ the one- loop RG analysis. 

By integrating out the fast modes, the one-loop RG 
equations are given by 
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(8) with the initial values: K a (0) = K a , K m (0) = K„ 



Do = K n 
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Ai(3)(0) - 7r fllC3) /(2« ), A 2 (0) = 9 2/(2v ), and A 4 (0) = 0. 
By solving Eqs. (|l(Jf) — l|15|). one may find three kinds 
of behaviors of the RG flow of Ai(Z) and X 3 (l): (i) Xi(l) 
flows to strong coupling while X 3 (l) flows to zero. This 
is the Id ASF. (ii) X 3 (l) flows to strong coupling while 
Ai(0 flows to zero. This is the I-CDW. (iii) Both Ai(Z) 
and X 3 (l) flow to zero. This is the 2-LL. Thus, Region IV 
in Fig. 13 (a) shrinks to a transition line between the Id 
ASF and I-CDW in the K a -K m space as shown in Fig. [3 
(b). Both the phase transition between the Id ASF and 
2 LL and that between the I-CDW and 2-LL belong to 
the KT type. The transition between the Id ASF and 
I-CDW is of second order. Further, all these transition 
lines coincide at two tricritical points, the point A and 
B in Fig. [21 (b). The very reason why there is no way 
for both Ai(/) and X 3 (l) flowing to strong coupling si- 
multaneously is that the operators cos ypHiQm — 2# a ) and 
cos v / 47r(</ ) a — 0m) are exclusive to one another, that is, 
the field configurations which minimize one perturbation 
term do not minimize the other. The interplay between 
these two competing relevant operators then produces a 
novel quantum phase transition. 

In the I-CDW, it is convenient to define new bosonic 
fields: (j) + = \{(j> a +<t>m), 4>- = <i>a-<i>ra, 0+ = a +6 m , and 
0- = \{0 a — 9 m ). The (f> + and 0_ fields describe the in- 
phase and out-of-phase density fluctuations, respectively. 
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Due to the relevant perturbation cos V47T(/>_ , the value 
of (4>-) is pinned and a gap is opened for the <fi- sector, 
while the 0+ sector is still gapless. On account of this, 
both the single-particle Green functions of atoms and 
molecules decay exponentially. On the other hand, the 
^Pa/b parts of the density fluctuations for atoms and 
molecules are enhanced: 



K m 
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< 2 for 7r|ff 1 |/(2« ),|fl 2 |/(2« ) 
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To understand the nature of the ground state of the 
I-CDW, a simple picture can be obtained in the limit of 
strong atom-molecule interactions (53), where the poten- 
tial energy (the g 3 term) dominates over quantum fluctu- 
ations. In this case, atoms and molecules form a regular 
lattice (Wigncr crystal of hard-core bosons). For g 3 > 0, 
the energy of the repulsion between atoms and molecules 



g 3 cos (V^ncf)-) = -g 3 cos [V47r(c/> a - <f> m + Vtt/2)] , 

is minimized by a relative phase shift of \pK jl between 
atoms and molecules, which corresponds to a shift of the 
atom (or molecule) lattice by half-a-period. 9] Thus, the 
I-CDW respects the symmetry of translation by one site, 
at — > di+i and bi — > 6j+i, but spontaneously breaks the 
reflection symmetry about the origin, at —* a_,; and 6,; — > 
b-i (or <j) a/m -> -4> a / m and 9 a/m -> 6 a/m ). 

Experimental signatures We suggest a few methods 
to detect the above results experimentally. The most im- 
portant distinction between the 2-LL and the Id ASF 
is the behaviors of the density correlators. For the 2- 
LLs, the 27T/9 a ( b ) part of the density correlators of atoms 
(molecules) exhibits power-law decay. However, the cor- 
responding sector decays exponentially for the Id ASF. 
Therefore, a threshold behavior will be observed at the 
momentum k = 2np a ^ in a time-of-flight measurement 
for atoms (molecules) in the Id ASF, while such a be- 
havior does not exist in the 2-LL. The other unique fea- 
ture of the ASF is that the ground state is a coherent 
state formed by hybridizing the atoms and molecules. 
Therefore, by suddenly changing the detuning, some kind 
of Rabi oscillation will be observed between the atomic 
and molecular condensates. Such a change of detun- 
ing can be achieved by applying magnetic field pulses 
to the ASF.0] Although there is no true condensate in 
Id due to strong phase fluctuations, a similar oscillation 
can also be observed between the densities of atoms and 
molecules, which is strongly damped by the phase fluctu- 
ations. Further support to the Id ASF is the examination 
of Eq. |(HJ), which can be achieved by a measurement of 
single-particle Green functions of atoms and molecules 
through the absorption line shape. 

Discussions Finally, two points should be addressed 
here. First of all, when ao\S\ is close to some rational 
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FIG. 3: Phase diagram for ao|5| being some rational number. 



number, one must replace the g 3 term in Eq. (|2J) by the 
g 3 term [Eq. J3J|], with the scaling dimension D2L A 
tree-level scaling analysis suggests the phase diagram in 
the K a -K m space as shown in Fig. |21 We note that 
there is no direct transition between the Id ASF and I- 
CDW. Moreover, the I-CDW can be reached only with 
very small values of K a i m . Next, a RG study shows 
that small velocity anisotropy, i.e. v a ^ v m , just shifts 
the phase boundary, and does not affect our conclusions. 
This is because the OPE's between the additional op- 
erators arising from velocity anisotropy and the already 
existing operators do not generate themselves. 

When this work was finished, we were aware of two re- 
cent papers dealing with a relate d sy stem — a Fcshbach- 
resonant atomic Fermi gas in ld.lllj Their results about 
the charge sector are similar to ours for the Id ASF. 
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